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Abstract 



We argue that there is a potential flaw in the standard treatment of weak 
decay amplitudes, including that of e'/e. We show that (contrary to conven- 
tional wisdom) dimension-eight operators do contribute to weak amplitudes, 
at order Gfo. s and without 1/Myy suppression. We demonstrate the existence 
of these operators through the use of a simple weak hamiltonian. Their con- 
tribution appears in different places depending on which scheme is adopted 
in performing the OPE. If one performs a complete separation of short and 
long distance physics within a cutoff scheme, dimension-eight operators oc- 
cur in the weak hamiltonian at order Gpo s //i 2 , fj, being the separating scale. 
However, in an MS renormalization scheme for the OPE the dimension-eight 
operators do not appear explicitly in the hamiltonian at order Gfol s . In this 
case, matrix elements must include physics above the scale \x, and it is here 
that dimension eight effects enter. The use of a cutoff scheme (especially quark 
model methods) for the calculation of the matrix elements of dimension-six 
operators is inconsistent with MS unless there is careful matching including 
dimension-eight operators. The contribution of dimension-eight operators can 
be minimized by working at large enough values of the scale /i. We find from 
sum rule methods that the contribution of dimension-eight operators to the 
dimension-six operator Qy is at the 100% level for fi = 1.5 GeV. This sug- 
gests that presently available values of \x are too low to justify the neglect of 
these effects. Finally, we display the dimension-eight operators which appear 
within the Standard Model at one loop. 
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I. INTRODUCTION 



The starting point for the study of nonleptonic weak transitions is the analysis of short- 
distance effects using perturbative QCD. The results are expressed using the Operator Prod- 
uct Expansion (OPE) as a series of local operators. In practice, only operators of dimension 
six are considered. We will argue that operators of dimension eight are also relevant, and 
that most previous analyses of nonleptonic amplitudes must be reconsidered. 

Nonleptonic weak amplitudes represent probably the most difficult calculations in QCD. 
Since the W-boson propagator is a constant up to Q ~ M w in momentum space , the 
amplitude is sensitive to strong interaction physics at all energy scales. Therefore, one must 
control simultaneously the very low, intermediate and high energy portions of the calculation. 
There are two key ideas, both introduced by Wilson [|T]-rj] , that are used in this regard. One 
consists of separating the different energy scales and integrating out those effects from high 
energy. This yields an effective low energy theory with modified interactions. The second 
is the tool for doing this - the Operator Product Expansion (OPE) - in which the effects 
of high energy are replaced by local operators ordered according to increasing dimension. 
The latter can be applied at any scale, and as we reduce this scale we successively integrate 
out more and more physics, thereby changing the coefficients of the operators. Specifically, 
one can consider the physics above and below some energy scale /x. Throughout the paper, 
we will refer to /i as the separation scale. For example, if one takes the complete set of all 
local dimension-<i operators^ {Q^} with the right quantum numbers, the operator product 
expansion tells us that a AS = 1 amplitude can be written to leading order in dimension as 

(n^ 1] ) = % WEE cfV) ■ (i) 

Here the {Q (//)} are coefficients which describe the short distance physics with Q > //. 
The subscript '/i' on the operator matrix element indicates that the matrix element is to 
include all physics up to the energy scale n {i.e. with Q < jj). The short distance OPE does 
not by itself solve the problem of nonleptonic amplitudes. However, it does tell us that the 
remaining task is to calculate the low energy matrix elements of local operators. 

An example of this appears in Fig. [I], where the long-distance (low-energy) and short- 
distance (high-energy) parts of a nonleptonic weak transition are depicted separately. Let 
the separation scale be /i. Then the gluons shown in Fig. |l|(a) have Q < \i and are associated 
with long-distance propagation. The blackened disc denotes all the (short- distance) effects 
with Q > fi. We next look into the short- distance regime via Fig. |l](b). Now there are only 
hard gluons with Q > fx which propagate over short-distances. If the separation scale /i is 
large enough, the physics at the higher energies can be analyzed with perturbative QCD. 
Because of this, the short-distance effects will appear local when viewed by low energy 
probes. Also, observe in Fig. [T](b) that the W-boson mass has been taken to infinity and so 
the process shown there corresponds to the range \i < Q <C M w . 



Throughout this paper, we explicitly display the operator dimension as a superscript, e.g. 0^ d \ 
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(a) (b) 
FIG. 1. Scales in the weak transitions: (a) Long range, (b) Short range. 

The most obvious candidates for the basis of weak nonleptonic operators are those with 
dimension six, formed as the product of two currents, 

q 1 Tq 2 g 3 rg 4 • (2) 

Since the hamiltonian is the product of two weak currents, the coefficients of these operators 
are dimensionless. Next will come the operators of dimension eight, some examples of which 
are 

QiFDpqa qaTWq* , (3) 
where is the covariant derivative and 

/ote?iyr^ 2 g 3 yr-g 4 i^ , (4) 

where F° is the gluon field-strength tensor. The coefficients of these operators have engi- 
neering dimension (Energy) -2 . In the original papers |5|,[J on the weak interaction OPE it 
was stated explicitly that the coefficients of dimension-eight operators are of order 1/M^, 
and this has been accepted ever since. All the current treatments consider only operators 
of dimension six 0-0. We will show that the correct procedure is more subtle. 

Specifically, we will first consider the situation where one has a true separation of scales 
in the fashion outlined above. In this case: 

1. Dimension-eight operators enter the weak OPE at order l/fi 2 and not 1/Myy. 

2. In the one calculable example that we know about, such effects continue to be signifi- 
cant even above the scale \i ~ 2 GeV. Since most present calculations are performed 
with fj, ~ 0.7 — > 2 GeV, dimension-eight effects will likely affect the results of past 
work. 

3. Most generally, the dimension-eight effect can appear in both the coefficient functions 
and the matrix elements. The relative amount of each depends on how one implements 
the division of physics at the scale \x and amounts to a 'separation scheme' dependence. 

However, in the process of demonstrating these points, we will also see that dimensional 
regularization does not accomplish the separation of physics above and below the scale \i. 
Since dimensional regularization is by far the easiest calculational scheme, we study the 
structure of the OPE in such a scheme. We find: 
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4. 



Matrix element evaluations in dimensional regularization must be sensitive to energies 
above the scale \i. 



5. In this case, the effects of dimension-eight operators appear fully within the matrix 
elements of dimension-six operators. 

6. Mixed evaluations, in which one calculates the coefficients using dimensional regular- 
ization and the matrix elements using a form of a cutoff, are inherently inconsistent. 
Most past calculations fall in this category. 

7. The influence of dimension-eight effects can be controlled by working at sufficiently 
large \i. Further work will be required to understand just how large /i must be to 
achieve a given precision. 

This work has two basic parts. In the first, we use an explicit analytic calculation to 
illustrate the properties of dimension-eight operator effects in a weak amplitude. This will 
provide a demonstration of the above points. In the second part, we calculate the relevant 
dimension-eight operators for the Standard Model AS = 1 weak hamiltonian in a particular 
separation scheme. This will allow the exploration of the size of such effects, provided the 
operator matrix elements can be evaluated on the lattice. 

II. AN EXPLICIT EXAMPLE 

Rather than deal with the usual weak hamiltonian, we start with a similar but distinct 
operator that has simplified properties and allows us to demonstrate analytically the exis- 
tence and properties of the dimension-eight operators. This hamiltonian [p!0| , pj]| contains one 
left-handed and one right-handed current instead of the usual Standard Model hamiltonian 
in which both currents are left-handed. Specifically we define^ 



H LR = ^J d'x 2V(x,M^) J^(x) , 
d(x)Tlu(x) u(0)T u R s(0) 

(Vf_ i2 (x) + ASLoCx)) (V7 +J5 (0) - A 4 ^ 5 (0))1 , (5) 



J^(x) = l -T 



= -T 

2 

where is the W^-boson propagator and V^, A% (a = 1, ... 8) are the flavor-octet vector, 
axialvector currents. 

The reason why this hamiltonian provides a useful example is that in the chiral limit 
its matrix elements are related to vacuum matrix elements, and we may therefore take 



advantage of what is known about the associated vacuum polarization functions. ||T2|,|T3j For 
example, the K-to-pi matrix element 



2 We omit CKM dependence in the operator Wlr and define the chiral matrices T^ L = 7^(1 ±75) 
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M(p) = (n-(p)\H LR \K-(p)) 



(6) 



is given in the chiral limit of zero momentum and vanishing light-quark masses by the 
vacuum matrix element 



<fx V(x, Ml-) (0\T (Vf(x)V^(O) - A» 3 (x)A r M) |0> . (7) 



Two of us have recently studied the amplitude M. and more details on its properties can be 
found in Ref. [H|. Here we will display those features useful for understanding the role of 
dimension-eight operators. 

One can perform an operator product expansion on the hamiltonian 7Ylr i n the usual 
fashion. Including only the dimension-six operators one findsQ 

— \cf\^)» + cf\^0P)»] . 



M 



2V2F? 



(8) 



The operator basis consists of two left-right operators Of\ which have respectively 
color-singlet and color-octet structure, 

T 3 . =. a. '■■'< 



(6) 



(9) 



In the above, q — u,d, s, r 3 is a Pauli (flavor) matrix, {A a } are the Gell Mann color matrices 
and the subscripts on 0[ , O^ refer to the color carried by their currents. The coefficient 
functions, including renormalization group summation, are 

8/9 / / \ \ -V 9 " 

I + 



1 



4%) = i 



a s (M ] 



w 



with 



a s (M w )) 



8/9 



\a s (M w ) 
a s {M } 



-1/9 



w j 



1 + 9- 



47T 



a s (M ] 



w i 



For our purpose it is sufficient to work with an expansion of Eq. 
a s {/i), 



M 



G f 



2V2F* 



In 



57T 



(10) 
(11) 

through first order in 
(12) 



Our goal is next to carry out an explicit evaluation of M. and to demonstrate that dimension- 
eight operators appear in addition to those of dimension six. 



3 We stress that {O k } and {of, } of this section are distinct from {Q\ } and {C\ } of Eq. (|l]). 
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A. The Presence of Dimension-eight Operators 



Let us analyze the vacuum matrix element that appears in Eq. (|7|) in terms of the vacuum 
polarization function 



(13) 



d'x e**{0\T(yf{x)V?(p) - A%(x)A»(0)) |0> 
Using this we transform the spatial integral in Eq. ((7|) to momentum space, 



M 



3G F M ] 



2 

IV 



32v^tt 2 F| 



dQ 2 



Q 2 + M^ 



U V ^Q 2 )-U A<3 (Q 2 



(14) 



Next we implement the separation of scales for the operator product expansion. We do 
this by applying a cutoff at Q 2 = fi 2 and using the OPE in the high-energy /short- distance 
portion. 

Thus consider a partition of M. characterized by the scale /i, 



M = M < (fi) + M>(fj) 



(15) 



where A^<(//) and J\4 > (n) are dependent respectively on contributions with Q < ji and 
Q > fi. We obtain then for the low energy portion, 



-m<(aO 



32V2tt 2 F2 Jo 
3Gp 



dQ 2 



Q 2 + M 2 V 



n y , 3 (Q 2 )-n Ai3 (g 2 



dQ 2 q 4 n y , 3 (Q 2 ) - n Ai3 (Q 2 ) + o(fi 2 /M } 



(16) 



32\/27r 2 F2 Jo 

This cut-off is well-defined as it refers to the external momentum of a gauge-invariant am- 
plitude. We have shown in Ref. [1C] that this relation serves as a definition of the vacuum 
matrix element for the local operator at the scale jj with a momentum-cutoff scheme 
(denoted as '(c.o.)'), 



dQ 2 Q 4 n y , 3 (Q 2 ) - U At3 (Q 2 



(17) 



More interesting for our purposes here is the high energy portion. We require that fi lie in 
the pQCD domain, and we further constrain it to obey fi <C M^. The asymptotic behavior 
of the vacuum polarization operator is then described by the operator product expansion, 
involving a series of local operators ordered by increasing dimension. In the chiral limit the 
leading contribution to the difference of vector and axial-vector correlators is a four-quark 
operator of dimension six |JI~2|fL4]| , followed by a series of higher dimensional operators, 



;n y , 3 -n A3 )(Q 2 



+ 



Q 6 Q 8 

Here £^ represents the combination of local operators carrying dimension eight. 



(18) 



These have been discussed and partially calculated by Broadhurst and Generalis [15 



For our purposes, it is not necessary to know their specific form, but only the fact of their 
existence. Upon performing the integration over Q 2 at high energies, we find 



M>{n) 



3G, 



In 



27r(a a Of) M +^ + 



P(8) 



In this expression we have dropped corrections of order fi 2 /M ] 
The full amplitude is then 



(19) 



M 



Gi 



2V2F 2 




(20) 



The crucial features here are the presence of the dimension-eight operators in the short 
distance portion of the amplitude M. and the fact that they appear divided by the scale 
fi 2 instead of M^. They are not suppressed by inverse powers of M w because these oper- 
ators appear in the vacuum polarization function at any Q 2 between fi 2 and oo. From this 
calculation it is clear that these operators must be present in an OPE that describes the 
integrating-out of short distance physics. These operators have previously been missed in the 
usual treatment of the operator product expansion within the weak hamiltonian, although 
they have been properly included in the OPE for the vacuum polarization functions |12[ . It 



will be clear from the work that we do below, where we find dimension-eight operators in 
the weak OPE, that they appear whenever we separate physics above and below the scale \i. 
Finally, despite our emphasis on the dimension-eight operators in this paper, operators of 
even higher dimension could play a role for sufficiently small values of \x. For example, the 
next term in Eq. (f^j) would be 3^ 10 V(327rV). 

where £^°> represents the dimension-ten 

effect. 



B. Estimated Size 

Here we give some numerical estimates of the size of the dimension-eight effects. This can 
be done using experimental data since the vacuum polarization functions satisfy dispersion 
relations. The inputs to the dispersion integrals, i.e. the imaginary parts are known from 
experimental work on cross section measurements of e + e~ — *> hadrons and from the study of 
hadronic final states appearing in r decay. We have performed the required phenomenology 
in Ref. ]nj, and the reader is referred to that work for more detail. Here we use that 
reference to illustrate the size of various effects in the OPE. 

First consider the sizes of the asymptotic elements in the vacuum polarization functions. 
Referring back to Eq. (0), we can display the relative size of the coefficients of Q~ 6 and 
Q- 8 . We find this to be 



£ 



(8) 

M-2 GeV _ i k 



-1.5 GeV' . (21) 



2n(a s o£ ) )^ 2 GeV 

Thus the dimension-eight effect is quite relevant for the \l ~ 1 — > 2 GeV region. 



6 



Let us also look at the magnitude of the three terms in the OPE shown in Eq. (|20|) . In 
the same order as displayed there (and in units of 10~ 7 ) we find 

r -0.12 - 3.84 + 0.64 + . . . (/i = 1 GeV) 

! n r rpV -2 M -0.28 - 3.49 + 0.30 + . . . fa = 1.5 GeV) 

iu i,ev yw - < _ Q 44 _ 3 24 + 1?+ (^ = 2 GeV) [ZZ) 

[-0.89 - 2.63 + 0.04 + ... fa = 4 GeV) . 

We see that at /i = 1 GeV, the dimension-eight term is larger than the leading operator 
Of ] in the OPE. fj Even at fi = 2 GeV, it remains a significant size relative to this operator. 
However, at \i = 4 GeV it is clearly small enough to be neglected. 

Through a great deal of effort, the short- distance perturbative structure of the weak 
interactions has been studied through two- loop order 0,§]. The dimension-eight effects are 
large enough that their neglect would negate this effort, as we would be left with only crude 
evaluations, at least at values of \x which are presently used. 



C. Separation Scheme Dependence 

If one changes the separation scale, there is mixing between the operators of dimension 
six and dimension eight. For example, if we use the scale \i\ instead of the scale fi, the local 
operator at the new scale is related (in our perturbative treatment) to those at the old one 
by 

(o! « )MHof » ) „ + _L ln g) (as0 ., > + g(J J _J ? ) . (23) 

Thus portions of the dimension-eight effect will appear within the dimension-six operator 
evaluated at a given scale. 

There is also a dependence on the scheme by which one performs the separation of scales. 
In the above example, we used a sharp cutoff in the variable Q 2 as the method of dividing 
the low and high energy regions. This is certainly the most convenient method in the context 
of the present calculation, yet need not be the only possible method. Imagine a smoother 
cutoff F(Q 2 //i 2 ), with F(Q 2 //i 2 ) -> 1 for Q 2 < /i 2 and F(Q 2 //i 2 ) -> for Q 2 > /i 2 . We 
assume that this function is such that all the following integrals are well behaved. Then let 
us define the vacuum matrix element of in a so-called 'F-scheme' as 

(o[x ] - jt~ 2 r dQ2 qa f (f ) [ ny - 3(g2) - nAs(g2) i (24) 

We also define the integrals 



*We will see that the appropriate comparison is with effect of oj 6 ^ rather than . 
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cp = — In 

dp 
7 



'3l 



+ 



M 2 W 



dQ 2 

g 2 q 2 + m 2 v 



l-F 



'91 



r°° dQ 2 


Hi?)] 


L g 4 





(25) 



We use these by inserting 1 = F + (1 — F) into the full matrix element, Eq. (|T4]). The first 
factor gives the of^ matrix element in the F scheme, and the remaining integrals can be 
done. In this scheme, the OPE for the amplitude A4 reads 

G p 



M 



2V2F% 



(6)\(F) 



(on 



3 
8^ 



In 



. ^ 2 . 



+ c F 



(oCsOP), 



+ 



3d} 



4 8) 



167T 2 /i 2 



+ 



Therefore the matrix elements in the two schemes are related by 



(6)\(c.o.) 



(on 



3cf 



(asO®), + 



3(l-d F )£® 



16tt 2 



+ 



(26) 



(27) 



To fully specify the OPE and the matrix elements, one needs to clearly specify the scheme 
for separating the scales. A lattice evaluation would involve different combinations than 
does our initial sharp cutoff scheme. We will return to this issue in the next section. 



D. Dimensional Regularization 

The presence of dimension-eight operators scaled by an inverse power of fi 2 appears odd 
in the method of dimensional regularization since one expects only logarithms of /i in that 
scheme. This is because in dimensional regularization one introduces an energy scale /id. r . 
('d.r.' denotes dimensional regularization) in order to maintain the proper dimensions away 
from d = 4. This scale appears only in the form /^~ 4 an d as d — > 4, fi^ r will appear only in 
logarithms. In this section we clarify this issue by means of explicit calculation. It will be 
seen that there exists a confusion between fi on the one hand and /id.r. on the other. The 
former is defined to be a separation scale (in the sense of effective field theory or of Wilson's 
OPE) whereas the latter has nothing to do with the separation of long and short distance. 
Indeed, dimensional regularization does not itself provide a mechanism for the separation 
of scales; all scales contribute to both the operator and the coefficient functions in such a 
scheme. 

We can evaluate the vacuum polarization functions as defined in d dimensions, 
l&i J d d x e^(0\T(Vf(x)V»(0) - A£(x)AZ(0)) |0> 

= W - fsTWvt - n A3 )(g 2 ) - gVn^ 2 ) . (28) 

^From this follows an expression for the vacuum matrix element of 0± in dimensional 
regularization, 

(0f)<£;> = (0|T (1^(0)^(0) - ^(0)^,3(0)) |0} 
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When d < 4, this expression is finite. A key point is that the integral continues to run over 
all Q 2 . It is not hard to relate this operator to the one found using a cutoff regularization. 
To this end, we split the Q 2 integral into regions below and above Q 2 = jj 2 . Note that 
fi is not the same as /id. r .- For the part of the integration below separation scale /i 2 the 
integral is finite for all dimensions, and we can take the limit d — > 4. This portion of the 
integration then reproduces exactly the cutoff version of the matrix element. We are left 
with the difference between the cutoff and dimensional regularization matrix elements. It 
comes entirely from the high energy region, i.e. with Q above the separation scale /i (see 
Ref. |TUJ for details; here we use the NDR scheme for 75), 



i-d 



(d - l)/i dj . 



dQ 2 Q d {u V)3 - n A3 ) (Q 2 



(47r) d / 2 r(rf/2) U 



(4n) d / 2 T(d/2) V 
3 



dQ 2 Q° 



P(8) 



Q 6 V ±) Q s 



+ 



57T 



A-d 



7 + ln(47r) + In 





r 







+ 



167T 2 yU 2 



+ . . . 



(30) 



In MS renormalization, the 2/(4 — d) — 7 + ln(47r) terms are removed, 
calculation one finds 



Completing the 



+ 



3as 
8tt 




P(8) 



167T 2 jJ? 



(31) 



Thus, all of the dimension-eight operator is shifted into the MS definition of the dimension- 
six operator. This is seen to be consistent: 

1. When one performs a separation of scales, one has the need for dimension-eight oper- 
ators in the OPE scaled by l//x 2 . 

2. When one defines instead the OPE using dimensional regularization, one cannot 
get effects proportional to l//x d _ r _, but the same effect appears contained within the 
dimension-six operator matrix element. Inspection of Eq. (BO) shows that one obtains 
the same overall matrix element^. 

It is important to recognize that dimensional regularization is not a true separation of 
scales for the OPE in the original sense meant by Wilson. |I|-|J The dimensionally regularized 
matrix element contains effects from all scales, including finite but sizeable contributions 
from short distances. Thus any use of dimensional regularization for the coefficient functions 
must be accompanied by an evaluation of the matrix element that covers all scales. 

The problem with this situation is that, in present practice, the matrix elements are 
always calculated with some form of a cutoff but the coefficient functions are calculated 



Since we are treating the dimension-six coefficients at leading-log order, we can ignore the non- 
logarithic dimension-six portion of Eq. (|3l|). To include it only requires an inclusion of the non- 
logarithmic terms in the coefficient function. 
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dimensionally. This is an inconsistent procedure. To relate operators in a cutoff scheme 
to dimensional ones, a dimension-eight effect needs to be included, as in Eq. (|3l"l) . We will 
return to a more complete discussion of this point later. 

III. DIMENSION-EIGHT OPERATORS IN THE STANDARD MODEL 

In this section we calculate the dimension-eight operators (and their Wilson coefficients) 
which are relevant for the Standard Model at one loop. We present these first employing 
a cutoff to provide a separation of scales, then return to a discussion of how to use these 
results in the context of dimensional regularization. 

A. Defining a Cutoff Procedure 

The construction of the OPE is performed by comparing a calculation performed both 
in the full Standard Model and within the effective theory. The coefficients of the effective 
theory are adjusted such that the results of the two are identical to a given order. For this 
purpose, calculations using free quarks and gluons are simplest, and are sufficient to identify 
the operators and their coefficients. 

We perform this calculation to one-loop order. For the full Standard Model amplitude, 
one calculates up to 0(a s ) in the usual way. For the effective theory one takes all matrix 
elements at tree level, except for those of the leading operator, 

Qf = uTjs dTtu , (32) 

which must be calculated at one- loop level in order to properly include the 0(a s ) effects. In 
both the full and effective theories one adopts the same external states and kinematics. It 
is important, however, to employ a distinct four-momentum for each of the external states. 

We use a cutoff to regularize the matrix element in the following way. Consider the 
matrix element of a current product at different spacetime points, 

M A ^ B = {B\T{jf{x)Jl^))\A) • (33) 

where J^ h is the hadronic charged weak current. Since the current is a color singlet, M.a-^b 
is invariant under QCD gauge transformations for any value of x. In order to define an 
operator matrix element from this, we fourier transform to momentum space, rotate to 
euclidean momentum and apply a cut-off on the euclidean momentum such that Q 2 < fi 2 , 

(B\0\A), = J |^ 9(/i 2 - Q 2 ) J d'x e^(B\T(J,(x)J^0))\A) . (34) 

This procedure puts the low-Q 2 components into the matrix element and leaves the high- 
Q 2 portion to be accounted for in the OPE. The analysis of the high-Q 2 portion is then 
accomplished in the same way as in QCD sum rules via the specification of the momentum 
flowing in the currents. In practice, this amounts to the following recipe. While calculating 
insertions in one-loop diagrams, we imagine the two currents to be connected by a VT^-like 
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boson with virtual four-momentum k. We route the virtual momenta in the loop according 
to this prescription, choosing k as the loop-integration variable and regularizing the theory 
by means of a sharp cutoff /i 2 on the euclidean squared- momentum k\. The result obtained 
this way is UV finite, and we define it as the matrix element at the scale \x. 

Given this separation of scales, the OPE operators can be calculated in one of two 
equivalent ways. 

1. We first calculate a matrix element in the full theory, which will be finite and inde- 
pendent of fi. The low energy radiative corrections to the operator below the 
scale n are then calculated in the effective theory While the infrared portions of these 
matrix elements will be the same, a comparison of the two reveals that specific local 
operators need to be added to the effective theory in order to reproduce the results of 
the full theory. These are the new operators in the OPE at one loop. In this method, 
dimension-eight operators are needed because radiative corrections of the matrix ele- 
ments in the effective theory contain l//i 2 effects, which do not occur in the full theory, 
and hence must be corrected for by dimension-eight operators. 

2. The same result can be found in the high energy portion of the calculation. In this 
method, the portion of the full theory that occurs above the scale fi is considered. 
This portion is equivalent to a set of local operators and the calculation readily reveals 
their coefficients. ]16[ In this case, it is seen that dimension-eight operators are a real 



contribution to the matrix element, and the factor of l//i 2 arises simply as the lower 
end of the region of momentum being considered. It is this method which was used in 
the first portion of this paper. 

These two methods are equivalent because the net physics is independent of \i. 



B. Results 



We summarize our calculation of dimension-eight contributions to both the current- 
current operators and the QCD penguin vertex. Throughout we work in the chiral limit 
of m u = rrid = m s = 0. This is an appropriate and useful limit because it captures the 
leading chiral contribution to kaon matrix elements. The leading weak chiral lagrangian 
that contains factors of the quark masses can be diagonalized away by a chiral rotation. 
This implies that the effects of quark masses are suppressed by one order in the chiral 
expansion. 



Box Diagrams: First, we consider the current-current sector of the AS = 1 four-quark 
sector which arises from the box diagram. We depict in Fig. ^] both the 'full' and 'effective' 
descriptions for one of the four possible box-diagram contributions. Virtual quarks which 
appear within a box loop are given a (very small) common mass W to regularize infrared 
behavior. 

To zeroth order in QCD the current-current hamiltonian is expressible entirely in terms 
of the dimension-six operator Q,f (c/Eq. (|32|)), 
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(a) (6) (c) (d) 

FIG. 2. QCD corrections to the box: full theory (a)-(b), effective theory (c)-(d). 




(b) (c) 

FIG. 3. Gluon radiation from the box diagram. 



H 



cff 



= %V us V: d Qf . (35) 

No QCD V 2 



The inclusion of QCD to first order yields the familiar dimension-six component, which we 
express in the 'color-basis' as 



where 



and 



o,(6) _ G_f_ v T/ * 

rL (curr-curr) ^/^ us ud 



qP = dry ut; s (37) 



CM 6) +4 6 W , (36) 



Cf\») = l + 0{al), C^(/,) = -|- ln f#) • (38) 

2 7T \ fi z 



The set of dimension-eight current- current operators is constructed from four-quark prod- 
ucts, covariant derivatives and gluon field-strength tensors. The gluon field-strength tensors 
we shall employ in our analysis are 

F% = d^G a u - d v G^ - g-sfabcG^Gl , F^ v = -e^ UQ/3 F^ . (39) 

As such, there will be contributions from gluon emission graphs (see Fig. |3|). Also, in 
constructing the dimension-eight operator basis, we have used the quark equation of motion 
and current conservation in the chiral limit, 

1 -Vq = 0, V fl (q 1 T^q 2 ) = . (40) 
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The list of gauge invariant dimension-eight operators is thenQ 



o (8) 


— ti U 1 v h til a LL "T til v Up u U b til Q ti T - til ^6 


A Ti TV" 


I a u + ui v s 


AV v Ti Ti^m 


Qf 










Qf 


= r> dr^Wn + «r£p M s J 5* r> 








Qf 


= u vj5 v v v a s dry + uv u a v u v^s dry + uY u a s 




ry + uvis 


dT^V^u 


Qf 


= uv^ T» a s dv v r> + uv a v v^s dryyu 








Qf 
Qf 


= uV^ T u a s dT»V v u + uT»VpS d V v T> 
= g 3 5 ab F^> b [uT^s dT u u - uT^s dT a v u] , 






( 



(41) 

where Z> M denotes a left-acting operation and we define the convenient notation 

r^^r^^ 7M (i+ 75 ) . (42) 



We find the Wilson coefficients corresponding to the above dimension-eight operators to 
have the form 



(8) _ Os_ J_ (8) 
i ~ A 2 '* > 



(43) 



where the {rjf} coefficients are: 



(8) 

m 


(8) 
V2 


(8) 

m 


(8) 

Vi 


(8) 
W 


(8) 


(8) 

m 


5/3 


22/3 


8/3 


-1/3 


16/3 


14/3 


1/3 



(44) 



QCD Penguin Vertex: In the chiral limit, the general form for the QCD penguin effective 
vertex is 



^-(QCD-pgn) _ G F 



(8)^,(8) 



i=l 



(45) 



Here, we employ a two-generation approximation (no dependence on the virtual top quark 
in the penguin loop). The penguin vertex as it appears in the 'full' theory is displayed in 
Fig. f|. Observe the presence of self-energy graphs in Figs. U(b),(c). 

The d = 6 QCD-penguin operator and its Wilson coefficient are given (in the chiral limit) 

by 



6 It might appear that the list in Eq. ( f4l| ) is missing another possible dimension-eight operator, 
namely gsfabcF^^uT^s dT^u. However, one can show that in the chiral limt, this is expressible in 

terms of Qf and Q,f- 
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9 



(6) 



FIG. 4. QCD penguin vertex in the full theory. 




(d) 



o 



(6) 



C 



(6) 



93 



(4tt) 



4 , /i 2 + m 2 



— In 

3 



+ 2g c - 2^ 6 2 + ^ c 3 



(46) 



where is the QCD coupling constant and g c is the dimensionless quantity 

,2 



/x 2 + m 2 



(47) 



Proceding next to the dimension-eight component, there are two classes of local, gauge- 
invariant operators: with an s — ► d quark bilinear and one field-strength tensor (Op}) and 
with an s — > d quark bilinear and two field-strength tensors ({Op 8 /*} (i = 2,3)). The latter 
correspond to two-gluon emission as in Fig. f|(d). Thus we find 



O p 8 } = dT a u s V a V a V^ a 



e> (8) 



U P3 



id 
d 



2 ' 2 



(48) 



along with the d = 8 Wilson coefficients 

#3 1 



ego*) 

cgV) 

4 8 3 V) 



(47r) 2 /i 2 + m 2 
a., 1 



lm 2 2 4 2 2 , 1 4 
3^3^ + 3^3^ + 15^ 



3ir n 2 + m 2 

Q s 1 
37r /x 2 + m 2 



4 m 2 



10 



13 



- 3# c + — ^ - — £ + -s. 



3 // 

4 m 2 . , 

~-j+2g B --gZ--g, 
3 /i z 



6 



4 
3' 



1 

3 5c 



(49) 



C. Conversion to dimensional regularization 

As we have stated, the dimension-eight operators and Wilson coefficients of the previous 
section are defined in terms of a particular momentum cutoff procedure. It is essential to 
understand how these relate to those obtained via other possible calculational approaches, 
most notably dimensional regularization and MS renormalization. 

In dimensional regularization, there are no l//i 2 effects. Therefore all dimension-eight 
effects must be transfered from existing as explicit operators in the OPE to being contained 
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within some MS matrix element. This was previously illustrated in our sample calculation 
Sect. II as embodied in Eq. (|3TJ). In the present case we are looking at radiative corrections to 
the operator Q% and therefore at this order all the dimension-eight effects will be absorbed 
into the MS matrix element of . It has long been realized that to convert to MS from 
some form of cutoff in the evaluation of matrix elements, a mixing of dimension-six operators 
is needed. What is new in the present work is the realization that dimension-eight physics 
is also needed. In terms of the coefficients {C\ } calculated above, the relation is 

(g (6) ) ^ =(Q (6) ) e.o + Edj(Q (6) ) ^ + Ec (8) (0 (8) ) ^ 

i i 

Here di are mixing coefficients for dimension six which we do not calculate in this paper. 
However, the key new feature here is that the matrix elements of the dimension-eight opera- 
tors must be added directly to that of dimension six in order to form the MS matrix element. 
Other MS operators will also have similar relations involving dimension-eight effects when 
converting from cutoff schemes. We discuss this point in the Appendix. 



IV. CONCLUSION 

We have uncovered a basic problem with existing calculations of weak amplitudes. The 
coefficients in the OPE are calculated using dimensional regularization and hence do not 
include dimension-eight operators. However, all matrix elements are calculated with some 
variation of a cutoff and hence also do not contain the effects of dimension-eight operators. 
When one connects the matrix elements to those in the MS scheme, one needs to consider 
both dimension-six and dimension-eight operators as our example shows in Eq. ([JTJ). However 
in practice this is not done, and the result is then inconsistent. 

Lattice evaluations of operators [I7|| rely on the finite lattice spacing to remove all short 



distance physics. This is a true cutoff, as the effects of short distance physics is simply not 
present in the simulation. Most current lattice calculations identify the inverse lattice spacing 
with the scale /i that defines the matrix element, although there are alternative possibilities 
(see item 1 below for an example). Therefore, to convert from lattice regularization at any 
fixed lattice spacing to a dimensional scheme requires the addition of both dimension-six 
and dimension-eight operators. This can be done as part of the 'improvement' procedure 
pOfl which attempts to correct for the effects of lattice artifacts, including the effects of finite 
lattice spacing. On the lattice the effects of higher dimension operators can be more severe, 
as the lack of chiral invariance can generate a dimension seven operator [ST], whose effect 



also needs to be corrected for. The present state of the art does not yet include the effects 
of dimension eight operators, but in the future it should be possible to correct also also for 
these. 

Quark model (and also large N c ) matrix elements are in a far more difficult situation 



In these cases, low energy models are postulated, and the models only make sense below 
1 GeV or so. These models are treated with a cutoff of order 1 GeV, and cannot contain 
the required short- distance physics. While lattice methods are correctible to account for 
dimension-eight effects, the same appears doubtful in quark model methods. We then must 
conclude that there are very large intrinsic uncertainties associated with these methods. 
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There are several strategies for overcoming the basic inconsistency we have revealed: 



1. In principle, the best way is avoid the need for a cutoff in the matrix element cal- 
culation - to have a method which is sensitive to all scales. This can then directly 
yield MS matrix elements. One possible procedure is the lattice analysis of the kaon 



I?-parameter in Ref. [19], where the scale /Zd.r. is held fixed while the lattice spacing is 
varied. The various values of the lattice spacing are then used to perform an extrapo- 
lation to the continuoum limit, at fixed fid.v.- Although this was done to remove lattice 
artifacts, it also has the effect of removing the need for dimension-eight effects (since 
the separation scale in the matrix element is sent to infinity). A second procedure is 
the evaluation that we provided in the first half of this paper, based on sum rule tech- 



niques where one knows the vacuum polarization functions at all scales ||10| , |iTfl . In this 
method there is no obstacle to a full evaluation, because we can take the separation 
scale to infinity^. See also Ref. |22| for a calculation sensitive to all distance scales. 



In general, the lattice seems most suited for a systematic program of extrapolating 
matrix elements to very large separation scales. We caution, however, that this has 
not yet been done in all cases of interest. 

2. A less ambitious but still valuable strategy would be to work at values of separation 
scale /i large enough that the problems of higher dimensions are numerically insignifi- 
cant. Stated in a different way, at large enough [i the matrix element evaluation will 
already contain all the needed ingredients, and residual dimension-eight effects that 
are missing will be small enough. This option is likely only available for lattice regu- 
larization, where it is tied to the ability to decrease the lattice spacing. Quark model 
methods make sense only at low energy and cannot be extended to larger fi. The 
question in this case is how large /i must be for a result of a given accuracy. We have 
addressed this for our model hamiltionian above with the result that fi ~ 3 — 4 GeV 
appears to be required. If the lattice is used to give estimates of matrix element of 
the dimension eight operators calculated in Sect. Ill, these results could be used to 
provide a second estimate of the relative importance of these terms. An advantage of 
this method is that one does not need to be highly accurate or highly consistent. If 
the dimension-eight effects are small at a given /z, then we do not need to match them 
to the rest of the calculation. We instead use the estimated size to produce an error 
bar due to this effect. If the error bar is small enough, we need need not consider 
dimension-eight operators further. 

3. Alternatively, one might adopt a cut-off scheme for both the matrix elements and 
the OPE coefficients. This would bring dimension-eight operators into the OPE, and 
their matrix elements would also need to be computed. However, this option appears 
difficult to carry out to the accuracy that we desire. At one loop, a cut-off scheme 
is only moderately difficult. One needs to be careful to preserve gauge invariance 
and other symmetres. One also needs to find a scheme that is equally useful for 



This will be discussed more fully in a separate work by us. 
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the perturbative calculation of the coefficients and for the calculation of the matrix 
elements if one is going to match them consistently. This latter requirement is likely 
difficult. However, it will be extremely difficult to implement a cut-off scheme at two 
loops. The presence of a dimensionful parameter upsets the normal power counting, 
and requires great care in properly defining the operators beyond one loop. In addition, 
the separation of nested loops in a cutoff scheme is subtle and can potentially lead to 
troubles with gauge invariance and other symmetries. 

4. Yet another option would be to continue to follow the most standard practice (of 
using dimensional regularization for the Wilson coefficients and a cutoff procedure 
for the operator matrix elements), but to correct the matrix elements to include the 
short distance effects from dimension-eight operators. In this case one would calculate 
both dimension-six and dimension-eight operators in some cutoff scheme, then add the 
contributions together with the right coefficient to form the dimensionally regularized 
matrix element. This is what is done earlier in Eq. fl31~1). This allows one to use the 
extensive work that has been performed calculating the OPE coefficients. It remians 
to be seen if this procedure can be successfully carried out in all cases of interest. 

Much of the existing work in the field is done at low values of \x. Various quark model 
and large-N c methods use fi ~ 0.5 — > 1 GeV, at which scales these effects are apparently 
extremely important. These calculations must be considered to contain enormous uncer- 
tanties, at least until further work is done. Lattice calculations are typically carried out 
with jj, — 2 GeV, although there is recent progress at working at higher \i. At this scale, 
dimension-eight contributions appear to still be larger than other effects which are included, 
such as scheme dependence and two-loop evaluation of the coefficient functions. We look 
forward to future work that allows us to eliminate or reduce the undertainties that come 
from the presence of dimension-eight effects. 
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APPENDIX A: ON THE MATCHING OF CUTOFF AND DIMENSIONAL 
REGULARIZATION PROCEDURES 



In Sect. Ill, we utilized the dimension-six local operator Q2 in our discussion of higher 
dimension effects in the Standard Model. Of course, a similar analysis can be performed 
for other operators Q- 6 ^ appearing in the weak hamiltonian. Here, we consider briefly the 
issue of perturbative matching between cutoff and dimensional regularization approaches in 
a general situation. 

Let us first establish the notation. We let Ql (fi) represent the d = 6 operators nor- 
malized in any cutoff scheme at the scale [i and Qi^ MS {^ms) represent the corresponding 
dimension-six operator in the MS normalization. Further, we let denote a tree-level 
matrix element (5!,- = ( Q [ ) ^ trec ' ) ) and consider matrix elements taken between states with 
generic momentum p and quark mass m q . 

Then, calculating in a cutoff scheme, and keeping terms up to dimension eight will give 
generally 



(Qi ) 



(6)\(c.o.) 



5ij + a s ~jij In 



p 2 



S 



(8) 



(Al) 



Here 7^ is the usual dimension-six anomalous dimension matrix. The coefficients Cjj scale as 
and therefore are sensitive to the infrared (IR) behavior of the matrix element, 



— 2 —2 

p or m, 



while %j is associated with the behavior at the upper end of the integration domain. f?j 0, 
are finite terms depending on the specific scheme adopted for separating the scales. For 
finite jj, the matrix element (Qf^)^' '^ will itself be finite and could be used as the definition 
of the operator matrix element at scale /1. Of course, the operators defined in this way differ 
from the MS ones by a finite normalization. 

We now present the connection to the operators in the MS scheme. This requires an 
appropriate matching in which dimension-eight operators will appear. An example of this is 
provided by Eq. fl3l|) in the study of the LR hamiltonian. Consider then a matrix element 
of the d = 6 operator Q,f^ taken between four-quark states. Expressed schematically, the 
calculation as performed in dimensional regularization ('d.r.') gives 



/ (6)x(d.r.) 

//M.r. 



5ij + a s ^ij In 



/4 

p2 



Sj -\- Q s CijS 



(8) 



with 



- = 7 + ln(47r) 

e e 



(A2) 



(A3) 



The infrared coefficients are the same as in the case of the cutoff scheme. In this case, 
however, the terms proportional to 7^ disappear because in dimensional regularization the 
integration runs over all scales. Introducing an MS-subtracted operator, 



(6)\(MS) 



''j 



(g (6) ){ d.,) 



(A4) 



we have 
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'q(6)\(MS) 



5ij + a s ^ij In 



MS 



+ OC s f, 



d.r. 



5 



(6) 



ot s CijS ' j 



(8) 



(A5) 



The general form of the above mentioned connection at order a s is given by 



'q(6)\(MS) 



MS 



5ij - a s -yij In 



J 1 - 



i _ i fWa tc.o.\ 



MS- 



{Qf) { r ] -^(of) 



(A6) 



The right hand side of Eq. (^G) is constructed in such a way to be finite in the limit 
in which fi — > oo. The /y matrices contain finite parts which depend upon the particular 
scheme adopted both in dimensional regularization and in the cutoff regularization. In 
particular, f^ s depends on the scheme definition adopted for 75 away from dimension four. 
The coefficients 7^ govern the 'leakage' of dimension-eight operators into the dimension-six 
sector. They depend on the particular scheme adopted for separating scales and, as shown 
before, they can be obtained by a one-loop calculation in the cutoff scheme. 
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